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ABSTRACT 

The dynamics of relativistic magnetized mass loaded outflows carrying 
toroidal magnetic field is analyzed in the context of Pulsar Wind Nebulae 
(PWNs). Mass loading is very efficient in slowing down super-relativistic 
magnetized flows and weakening of relativistic shocks. We suggest that 
weakening of relativistic reverse shocks by mass loading in PWNs is responsible 
for the low radiative efficiencies of the majority of the PWNs. Mass loading may 
also result in a shock transition near the fast magnetosonic point; this is unlikely 
to happen in majority of PWNs. The evolution of magnetized mass loaded 
flows beyond the reverse shock is complicated: after initial deceleration to the 
minimal velocity required to transport the magnetic flux, the mass loaded flows 
have to accelerate. In order to be able to expand to infinity, magnetized flows 
should either become time dependent or destroy the toroidal magnetic flux by 
developing internal instabilities. Destruction of the magnetic flux initiated by 
mass loading may allow for the flow to slow down to sub-relativistic velocities 
and resolve the a paradox of the pulsar wind nebula. 



1. Introduction 



In this paper we explore the dynamics of magnetized relativistic pulsar winds subject 
to mass loading. We are motivated by several possible applications: (i) slowing down of the 
high-cx flows [a is conventionally defined as a ratio of Poynting to particle fluxes) inside the 
pressure confined pulsar wind nebulae (referred below to as a "<r paradox"), (ii) low X-ray 
efficiencies and a conspicuous lack of radio PWN around many energetic pulsars and (iii) 
the structure of the ram pressure confined pulsar wind nebulae. 

(i). The a paradox is a long standing problem in pulsar physics (Rees & Gunn 1974, 
Kennel & Coroniti 1984). Models of the pulsar magnetosphere (Goldreich & Julian 1969, 
Arons & Scharlemann 1979, Ruderman & Sathurland 1975) predict that near the light 
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cylinder most of the spin-down luminosity of a pulsar should be in a form of Poynting 
flux, a >> 1. On the other hand, modeling of the dynamics of the Crab nebula gives a low 
value of o at what is commonly believed to be a reverse shock - strongly magnetized flows 
cannot match boundary conditions (Kennel & Coroniti 1984). Several possibilities have 
been proposed to resolve the a paradox. One possibility is that o may change between the 
light cylinder and the reverse shock (Michel 1994, Coroniti 1993, Melatos & Merose 1996, 
critique by Lyubarsky & Kirk 2001). Alternatively, a may change between the reverse shock 
and the contact discontinuity. Begelman (1998) argued that the post-shock flow may be 
distorted by kink instabilities which would annihilate the magnetic flux and allow strongly 
magnetized flows to match the boundary conditions at the edge of the nebula (see also 
critique by Arons 1998). Arons (2001) has suggested that interaction of the wind with the 
ejecta filaments may slow down the wind. In this paper we investigate the latter possibility. 

(ii) . Most PWNs, with the notorious exception of the Crab, have very low efficiencies 
for the conversion of the pulsar spin-down luminosity Eq into X-rays, Lx ~ 10~ 3 Eq (Seward 
& Wang 1988, Becker & Truemper 1997). No simple combination of the pulsar parameters 
(like period and period derivative) can be fitted to the X-ray luminosity, which suggests 
that environmental effects may be important (Possenti et al. 2001). Similarly, sensitive 
radio searches have failed to detect PWNs towards selected energetic pulsars (Gaensler et 
al. 2000). We show that mass loading of the pulsar wind is very effective in slowing down 
the wind and thus weakening the termination shock. This may be responsible for low X-ray 
and radio efficiencies of the PWNs. 

(iii) In the past several years a new well defined class of objects - ram pressure 
confined pulsar wind nebulae - has been established. Apart from a simple model by 
Wang et al. (1993) no realistic model has been proposed to account for the often unusual 
structure structure of such PWNs. Bucciantini & Bandiera (2001) have argued that classical 
bow-shock models cannot be applied to the ram pressure confined PWNs and suggested 
that mass loading may be important. 

The importance of mass loading on the pulsar wind may be estimated using the 
Crab nebula as an example. It has long been known that that some thermal plasma does 
evaporate off the trapped ejecta (Wilson 1974, Michel et al. 1991). It is only natural to 
expect that the coupling of a "light" pair plasma by the "heavy" particles from the ejecta 
would strongly affects the dynamics of the light fluid. Recent reevaluation of the mass of 
Crab filaments (Fesen 1997) gives ~ 4.6 ± 1.8M , while the total mass ejected into the Crab 
nebula by the pulsar over its lifetime is tiny ~ 1O~ 12 M (at a rate N = 10 38 pairs/sec (e.g., 
Kennel & Coroniti 1984); ions may have a two order of magnitude larger mass flux; Galant 
et. al 1992). Below we show that mass loading starts to strongly affect the dynamics of 
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sub-sonic flows when the pick-up mass flux, M becomes of the order of the total energy 
flux E divided by c 2 (eq. flSS|)). The total pick-up mass flux is M ~ M picked v/r (with 
velocity given by the expansion velocity v ~ 2000 km/sec of the Crab nebula with the 
radius r ~ 2 pc; Kennel & Coroniti 1984). The mass flux required to slow down the flow to 
match the boundary M P i C k e d ~ Eqt/(c 2 v) = 3 x 10 27 g, a tiny fraction of the available mass. 
Thus even a weak coupling of filaments to the flow may strongly affect its evolution. This 
may be contrasted with the theoretical estimates of Balbus (1981) who have argued that 
filament evaporation by conduction may occur on a very short time scale, of the order of 
the dynamical time scale. 

In the ram pressure confined PWNs, we should substitute for v the velocity of the 
pulsar v ~ 200 km/sec and for r the stand-off distance, which is typically a fraction of a 
parsec. The total pick-up mass flux, M P i C k e d ~ 10 24 g, is nevertheless smaller since typical 
pulsars have Eq ~ 10 35 erg/sec. This is approximately the amount of the ISM material 
contained in the head part of the ram pressure confined PWN. 

2. Model assumptions 

We assume that an interaction between the trapped ejecta and the pulsar winds occurs 
through mass loading of the wind. Filament material may be coupled to the pulsar wind 
by thermal evaporation, hydrodynamic ablation and by photo- and electro-evaporation 
by high energy photons and relativistic electrons. Interaction of the ISM clumps with 
non-relativistic winds have been studied extensively (Borkowski et al. 1990, Cowie & 
McKee 1977, McKee & Cowie 1977, Hartquist et al. 1986). No treatment of the interaction 
of the ISM material with relativistic magnetized flow exist; estimates based on the models 
of suprathermal evaporation (Balbus & McKee 1982) predict a very high rate (Balbus 
1981). Here we are interested in the dynamical effects of the mass loading on the relativistic 
flow; hence we do not consider here the microphysics of the flow-filament interaction and 
instead we parameterize it by a capture rate, (eq. (§)), taken to be a prescribed function 
of a distance from a pulsar. The newly acquired particles are assumed to be immediately 
coupled to the flow. [] We neglect the energy losses required for the possible ionization and 
the energy associated with the thermal motion of donor particles. Mass loading inside a 
shock may also be important (Zank & Oughton 1991). Mass loading of relativistic shocks is 
briefly considered in the Appendix [A]. 



1 Relative streaming of the newly ionized component with respect to the bulk flow may result in ion- 
cyclotron instabilities and turbulence generation (e.g., Galeev et al. 1996) 
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The dynamics of the mass loaded flows turns out to be non trivial and somewhat 
counterintuitive. The overall dynamics of the mass loaded flows resembles a nozzle-type flow 
where the rate of mass-loading contributes a term which may be considered as a negative 
pressure (eq. (fZUD). Hence the effects of mass loading are opposite to that of pressure. 
For example, contrary to the naive guess, the loading of strongly magnetized pulsar wind 
leads not to slowing down but to acceleration of the flow. Acceleration formally proceeds 
to arbitrary larger velocities which, naturally, cannot be realized in reality. We argue that 
such acceleration of the wind will result in development of internal instabilities (which we 
have not addressed yet) that would destroy the toroidal magnetic flux and, in a way similar 
to suggestion by Begelman (1998), would eventually allow a wind to slow down. 

The pulsar wind model that we adopt is similar to the Kennel & Coroniti model: 
spherical outflow of a polytropic gas with a toroidal magnetic field. In addition, the wind is 
assumed to be loaded by ejecta particles which are initially at rest. By neglecting the radial 
component of the magnetic field we limit the applicability of this approach to asymptotically 
large distances, far from the acceleration region of the flow (Weber & Davis 1967, Goldreich 
& Julian 1970, Kennel et al. 1983). Typically the launching of the pulsar wind is thought 
to occur close to the light cylinder, with the flow reaching super fastmagnetosonic velocities 
by passing successfully through the slow, Alfven and fast sonic points (Goldreich & Julian 
1970). At distances large compared with the light cylinder radius, the flow usually can be 
well approximated as a radial motion with toroidal magnetic field, launched with some 
given parameters: mass, energy and magnetic fluxes loss rates. 

Pulsar winds are launched at super-fastmagnetosonic velocities and have to decelerate 
in a shock transition in order to be matched onto the non-relativisticly moving external 
media (the SNR ejecta or ISM). Under the hydrodynamic approximation, a contact 
discontinuity separates the shocked wind material from the external medium. In some 
circumstances, the external material may get inside the contact discontinuity and even 
inside the reverse shock propagating in the wind. This may be due to either incomplete ISM 
or ejecta ionization (neutral particles may not be well coupled to the ionized component 
on time scales of interest to us) or by evaporation of the ISM or ejecta clumps trapped by 
the expanding wind. Coupling of the newly added material to the flow may considerable 
change its evolution. 

The evolution of non-relativistic mass-loaded flows has been investigated in a number 
of papers, starting with the Solar wind-comet interaction (Biermann et al. 1967). Lately 
Hartquist et al. (1986), Arthur et al. (1993), Smith (1996), Williams et al. (1999), Toniazzo 
(2001) have investigated effects of mass loading on stellar outflows. Here we show that the 
presence of a magnetic field makes it formally impossible for the mass loaded flow to reach 
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infinity. In order to be able to reach infinity, the flow must either become non-stationary 
or destroy its magnetic flux, presumably by developing internal instabilities which result in 
reconnection. 



3. Governing Equations 

The formal treatment of the problem starts with the set of relativistic 
magnetohydrodynamic equations which can be written in terms of conservation laws 
(e.g. Landau & Lifshitz 1975): 

n = sk, (i) 

Ff = 0, (2) 

(jm'li = R (3) 

where 

T ij = (w + 6 2 )wV + (P + ^ 2 )g ij ~ V V (4) 

is the stress-energy tensor, w = p + rcrP is the plasma proper enthalpy, p is proper plasma 
density and p is pressure, b 2 = bJcP is the plasma proper magnetic energy density times 47r, p 
is pressure, u l = (7, j/3) are the plasma four- velocity, Lorentz-factor and three- velocity, g 1 ^ is 
the metric tensor, 6j = hrjiji-iU^ F kl are the four-vector of magnetic field, Levy-Chevita tensor 
and electro-magnetic field tensor, Si is a four-vector representing energy and momentum 
sources and R is the density source (dimensions of R are gem -3 sec -1 ). We assume that the 
adiabatic index T is constant for algebraic simplicity. Since we are interested in transonic 
transitions for relativistic fluid motion, we expect that a relativisticly hot flow with p^> p 
will have an adiabatic index of 4/3 while strongly subsonic (non-relativistic) flows will 
asymptote to T = 5/3. 

We assume that loading is due to the medium at rest in the laboratory frame, S = 0, 
S = R. 2 Writing out eqns. (Hl-^D in coordinate form and assuming a stationary, spherically 
symmetric outflow with toroidal magnetic field, we find 



r 2 {w + b 2 )(3^ 2 \ = R (5) 
r 2 ((«, + 6 2 )/3 2 7 2 + (p + b 2 /2))\ - ^ = (6) 



2 The filaments are moving with velocities up to 300 km/sec (Trimble 1968), much smaller than the 
velocity of the pulsar wind near the termination shock w s hock ~cor near the edge of the remnant v ~ 2000 
km/sec (Kennel & Coroniti 1984). 
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^d r [rbPi] = (7) 
' r 2 p/3 7 =R (8) 



The above relations can be simplified if one introduces energy, mass and magnetic flux loss 
rates, 

c = m 2 { b2 + Y~jP + p) ( 9 ) 

T = {3 1P (10) 
JC = [3-fb (11) 

which evolve according to equations 

2£ 

d r C = + R (12) 

r 

IT 

d r T = + R (13) 

r 

d r JC = -— (14) 



which have solutions 



r 



j~ _ Mp + M 
4 7r r 2 
_ E + M 
Anr 2 

£=7t4=- ( 15 ) 

where M and E'o are the central source's mass and energy loss rates, £ is the electromotive 
force and we have introduced 

M = 4vr J Rr 2 dr (16) 

for the acquired mass flux. 

It is convenient to introduce two other parameters of the flow: magnetization a and a 
fast magnetosonic wave phase velocity (5f and Lorentz factor 7/ 

b 2 K 2 £ 2 
a = — 



w Cp-K 2 E f3 - £ 



2 



l + o- (l + a)w 
, _ 1 3£(1 + a) 



7/ " 1-0} 2£ + ^ 7 (l + a) (17) 
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In particular for a relativisticly hot plasma, p <C p, we have w = Tp/ (T — 1) and 
^= " + r - 1 = 1 + 3(7 for T = 4/3. 

M/ 1 + 0" 3(1 + (X) ' 

We will need also an expression for pressure in terms of the fluxes: 

£-g7(l + g) _ (r-1) (£-^77/ 2 ) 

p 1 ' v^n + a) m 2 (2-r) r 7/ 2 

Since p > it follows that £ > 7 (1 + a) T and £ > 7 T . 



(19) 



Eliminating /C in favor of 7/ we get a particularly transparent form for the evolution of 
Lorentz factor 

( < 2 - r > ^ ( 7 3 - 7/ .) 8r7 = 2 < £ V) _ (7 - i) (2-r) (i + 2 r) V fl (20) 



\/P>f (r-i) j v " ' cr £(r-i) 

Equation (p0| ) has a form of nozzle- type hydrodynamical flows (e.g., Landau & Lifshitz 
1999); in astrophysical context it is best known for Parker's solutions of the solar wind 
(Parker 1960). The lhs of eq. ( p0|) contains a familiar special point at the sonic transition 
7 = 7/. The positively defined first term on the rhs describes the evolution of Lorentz 
factors due to pressure effects, and the negatively defined second term is due to mass 
loading. Thus, the rate of mass loading may be considered as a negative pressure. 

By neglecting the radial magnetic field far from the acceleration region we have lost 
the magnetic sling-shot effect often evoked for the acceleration of plasma (Michel 1969, 
Goldreich & Julian 1970, Kennel et. al 1983). The upshot of those works is that the flow 
is accelerated to supersonic velocities with a terminal Lorentz factor 7 ~ ^fo. Since it is 
believed that the plasma near the light cylinder has o"o ~ 10 3 — 10 6 , but the inferred pulsar 
wind Lorentz factor is ~ 10 6 3> y/&o, an additional acceleration is required. 



Without mass loading the rhs of eq. fl20|) is always positive so that super- 
fastmagnetosonic flows accelerate while sub-fastmagnetosonic flows decelerate. With mass 
loading the situation is more complicated. Mass loading enters in the eq. ( p0|) in two ways: 
explicitly via the rate of mass loading R, given by the negatively defined second term on the 
rhs, and implicitly via the accumulated mass and energy fluxes in T and C The negatively 
defined second term could in principle become larger than the pressure term. Below we 
show that this indeed happens for magnetized flows. 



A quick examination of eq. ( f20|) reveals an important fact: for strongly relativistic flows 
the mass loading term is enhanced by a factor ~ 7 2 . For typical pulsar winds 7 ~ 10 6 , thus 
mass loading is extremely efficient in slowing down the wind. We defer a detail examination 
of the eq (|20| ) until Section |5] and review beforehand the evolution of unloaded flows. 
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4. Dynamics of relativistic flows without mass loading 

To guide us through the effects of mass-loading we first derive the relations governing 
the evolution of unloaded relativistic magnetized winds. The properties of these solutions 
have been discussed extensively (e.g. Kennel & Coroniti 1984), but the exact analytical 
form given below has not been written, according to our knowledge. First, we introduce 
two new parameters (instead of M and £) - the unloaded wind terminal Lorentz factor 7 
and the initial magnetization parameter oq defined by the following reactions 

Eq 

7o 



M (l + (To) 

S 2 



a = (21) 

To/Wo 



4.1. Unmagnetized unloaded relativistic flows 

In the absence of a magnetic field the fast magnetosound velocity is 

0} = (T - 1) (£ ~/ 7) . (22) 

In this case the equation for the evolution of the wind Lorentz factor (|20|) can be integrated 
exactly: 

r oc _L ( 7o _ 7 )i/(2(r-D) 7 (2-r)/(2(r-i)) _ (23) 

VP 

This shows that there are two branches of solutions: supersonic and subsonic. The 
supersonic branch after initial acceleration with 7 ~ r 2 ^" 1 '^ 2 ^ (7 ~ r for T = 4/3) 
reaches a terminal Lorentz factor 70 = E /M (with 70 — 7 ~ r -2 ^ -1 )). In the accelerating 
part the flow has p ^> p - it is a pressure driven acceleration. Since pressure falls off with 
distance faster than density (for T > 1), a transition to the coasting phase with 7 ~ 70 
occurs at p ~ p. The subsonic branch decelerates to zero velocity at infinity with j3 ~ 1/r 2 , 
keeping pressure and density almost constant, determined by the mass conservation ratio. 



4.2. Relativistic Magnetized Unloaded Flows 

Magnetized flows also have two branches: subsonic and supersonic. The terminal 
velocity is determined from the condition d T r ) = = p, which, using eq. (|T7|) , can be written 

7-7o P-Po ^ (0A , 
= — n — ^o- (24) 

7o p 
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Eq. fl24]) generally has two solutions: a supersonic one 

7 = 7o (25) 
(it is highly supersonic for 7 ^> yfoo) and a subsonic one, which in the limit 70 ^Jo\ ^ 



gives 

P = T~rr- (26) 

1 + 00 

No physical solutions exist for (3 < 00/ (1 + 0q) - such flows cannot transport enough 
magnetic flux. 

The supersonic solution behaves similarly to the unmagnetized case. The flow is 
accelerated by pressure effects as long as p 3> p, reaching a coasting phase with 7 ~ 70 when 
p < p. Behavior of the subsonic branch is qualitatively different from the unmagnetized 
case. To study the behavior of the subsonic branch we can use a simplifying assumption 
7o — > 00 since for subsonic flows the Lorentz factor is usually only weakly relativistic, 
7 < \J (1 + cr ) /(2 — T) while 70 ^> y/cro- The limit 70 — > 00 is equivalent to neglecting the 
mass loss rate of the central source if compared with the energy loss rate. The evolution of 
the flow is then given by 

r oc (/^-^(-D) (/3 _ , o/(1 + ^-VCCr-D) = f 7 for L = 4/3. (27) 

(p - 0o/(l + 0o))' j/i 



This shows that subsonic flows reach a minimum velocity given by eq. (p6|). The fact that 
magnetized flows cannot slow down to zero velocity since they have to transport magnetic 
flux is crucial in determining the asymptotic dynamics of the mass loaded flows. 



5. Evolution of mass loaded pulsar winds 

A pulsar produces a wind with properties determined by three parameters: energy 
flux E , terminal Lorentz factor 70 and magnetization ctq. We also assume that near the 
acceleration region (located presumably close to the light cylinder) the mass loading is 
unimportant, so that the wind has reached a strongly super-fastmagnetosonic terminal 
Lorentz factor 70 3> \/l + cr (this requires M £ 2 /(?o, E (3 /(l + cr )) and its internal 
pressure has dropped to 0. 

The ratio of the mass flux due to loading to the initial energy flux is 

47r f Rr 2 dr AirRr 3 , . 

C = — —■ = 28 

E (3-n)E 

for a power-law dependence of mass loading on radius, R ~ r~ n . Parameter ( will often be 
used instead of the radius r. 
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5.1. Weakly loaded flow 

Here we show that mass loading of relativistic flows is extremely efficient. Consider 
evolution of a strongly super-fastmagnetosonic cold flow when mass loading may be 
approximated as a small perturbation. In the limit r — > and 7, 70 ^> 1, eq. (|20|) gives 



<91n 7 2(70-7) ^r 2 R r7 2 7 (l + a ) 



d\nr 7 + (2 - r) 7o /(r - 1) E "f(T - 1) + (2 - r) To 
If initially 7 = 70, then 



(29) 



9 In 7 47rr 3 i? 2 4nr 3 R 

ainT = -^T r7o(1 + ao) = "^r r7 ° (30) 

Thus, the flow starts to decelerate quickly when the accumulated mass flux becomes 

comparable to Mo/70 - an increase in efficiency of loading by a factor of 70 (~ 10 6 !). The 
typical scale for the deceleration of the flow is 

( E \ 1/3 / Mq \ 1/3 

Td ~ UnRT^l + *))) ~ WrT 1o ) {6 > 



1 /3 

which is 7 times smaller than in the non-relativistic case. 

Efficient deceleration of supersonic flows by mass loading come from the requirement 
that the newly acquired particles, which were initially at rest, have to be accelerated to 
large Lorentz factors. This costs a lot of energy and momentum. 



5.2. Special points in the flow 

Special point in the flow occur when the rhs of the eq. ( p0|) becomes 0. It is shown 
below that in unmagnetized flows, the only special point in the flow occurs where f3 = j3f, in 
order to keep d r /3 finite, while mass loading of magnetized flows introduces another special 
point where /3 ^ j3f and, thus, d r (3 = 00. 

The rhs of eq. ( p0| ) is equal to only at one point, where the lhs of the eq. ( pOf) is also 
0. This is a well studied critical sonic point of the flow where (3 = (3f. For a magnetized 
flow the rhs of eq. (|20| ) become equal to at two points: a critical sonic point and another 
point located for a given flow characteristic at radii larger than the sonic point. For even 
larger radii a sub-fastmagnetosonic mass-loaded flow starts to accelerate - this behavior is 
opposite to the unmagnetized flow. Acceleration of the flow proceeds very quickly; the flow 
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reaches a point f3 = (3f and d r (3 = oo. A steady state flow cannot exist beyond this limiting 
radius. These conclusion hold for both relativistic and non-relativistic flows (considered 
in the appendix [B] The phase portrait of the mass loaded relativistic magnetized flows is 
given in Fig. and non-relativistic in Fig. |^. 



5.3. Sonic point 



Sonic point occur when both side rhs of the eq. (|20|) are equal to zero. It is possible to 
find the location of the point (3 = (3f on the phase diagram (3 — (. The general relations are 
complicated, yet simple approximations may be obtained in the limiting case of small and 
larger o~o- I n the limit o~o <C 1 the sonic point is located at a Lorentz factor satisfying 

3 (2-r) r 7/ 3 + (4-r) 7/ 2 - (2 + r) 7/ -3 = o. (32) 

The position of the sonic point in radius is given by Cf = f(T) = 0(1), where f(T) is a 
complicated function of the order of unity. For example for relativistic flows, T = 4/3, 
7/ = 1.09, f3f = 0.40 and Cf = 0.91 while for non-relativistic flows, T = 5/3, 7/ = 1.31, 
/3 f = 0.65 and (/ = 0.39. 

In the opposite limit a ^> 1, we find 

2 4-r . , 

^ " W^n ao (33) 

and Cf = /(r) / a^ 2 where again / is a complicated function of the order of unity. This gives 

7 | = 2a , Cf = 1-19 x °o 5/2 for T = 4/3 

7/ = 5^0, Cf = 0-23 x a 5/2 for T = 5/3 (34) 

Thus, low a flow experiences a shock transition when the swept up mass flux becomes of the 
order of the luminosity of the central source, while high a flow experience shock transition 
at a distance ~ crj^ 6 closer to the source. Location of the sonic points for two choices of do 
are shown in Fig. 0. 

Since in the hyper-sonic regime 70 falls out of the equations, flows with different 7 o but 
the same o"o experience a shock transition at the same location in C coordinate (Fig. IT]). 



5.4. Heavily loaded unmagnetized winds 



For unmagnetized flows the effects of mass-loading may completely dominate the 
dynamics of the flow at larger distances. This occurs when the loaded mass flux becomes 
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larger than the initial mass flux. In this limit M 3> E Q , M Q and £ = it follows from the eq. 
( pOD that only non-relativistic strongly loaded flows can extend to infinity Asymptotically 

„ r -(6T+l)/(r+l) ^ r -7/2 (35) 

for T = 5/3. 



5.5. Evolution of the magnetization parameter a 

Finally we consider evolution of the magnetization parameter a. Restricting to T = 4/3 
we find 

7 2 ~7| a 1-.T77? 7(7-1 + 2/?^ D 

n , \ dr<? = ? + R 36 

a(l + a) Lr 2£ 

with Pf and 7/ given by eq. (|T7|). Thus effects of mass loading are opposite to that of 
pressure. In the absence of mass loading, pressure effects lead to a a decrease for supersonic 
flows and a a increase for subsonic flows. Mass loading contributes to an increase in a for 
supersonic flows and decrease for subsonic flows. 

In mass loaded supersonic flows, the pressure term always dominates over the mass 
loading term in the eq. (|3T)|), so a always increases, formally diverging at the point (3 = f3f. 



6. Discussion 

In this paper we have investigated the dynamics of relativistic magnetized mass loaded 
flows. We found that (i) super-relativistic flows are effectively slowed down by mass loading; 
(ii) subsonic magnetized flows subject to mass loading have to accelerate beyond some 
radius and formally cannot reach infinity. We suggest that internal instabilities would 
develop in the accelerating flow, destroying the toroidal magnetic flux. Alternatively, the 
flow may become non-stationary. 

To estimate whether mass loading is important in slowing down pulsar winds we need 
to know the loading rate R which depends on the microphysics of the flow-ejecta coupling. 
To make a simple estimates we assume that a considerable fraction of the ISM neutrals 
(or trapped ejecta particles) with density (1 — £) n o> ionization fraction £ and mass m p are 
captured by the flow on a characteristic time scale r: R ~ ra m p /r. Since in the Crab we do 
see filaments (i.e., they are not evaporated on a dynamical time scale), for static PWNs the 
age of a pulsar T may be considered as a lower limit onr >T and thus an upper estimate 



of the rate R. From the eq. (31) it follows that the flow's initial Lorentz factor decreases 
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on a scale 

r „ 10 « (J^S" (^4) m (-Z-f f^r" fil-)- 1 ' 3 c,„ (37) 
\^10 36 ergy V 0.1cm" 3 J \10 4 yrs/ UOV V 1001 / v ; 

We would like to stress that this is a scale for a change of the initial Lorentz factor and 
not the scale for the mass-loading-induced shock transition, which is much larger (see eq. 
(|38|). Distance fl3~T| ) has a weak dependence on the characteristic time scale r, so that the 
deceleration of relativistic winds due to pick-up ions may be an important factor affecting 
the evolution of pulsar winds. 

To estimate the position where a mass-loaded shock would occur we assume that 
pick-up rate is independent of radius. Then the shock transition occurs at 



3E t \ 1/3 1 , for (T < 1 



r ~ — — x K/a ' ~ 10 x FTTrcm (3 

\^{l-i)n m p ) a 5/6 , for a > 1 1 + J/* V 



(for the same parameters as in eq. (pTf)). Radius (qQ) also gives a typical location where a 
flow becomes dominated by mass loading. Radius ([38|) is of the order of the shock distance 
in the Crab nebula. Since it is only a lower bound (since r ^> T) we conclude that a mass 
loaded shock transition is unlikely to occur in the static PWNs. The shock observed in 
the Crab nebula nebula is likely an ordinary magnetohydrodynamical shock. On the other 
hand radius (p8|) is much smaller than the radius of the Crab nebula. Thus mass loading of 
the subsonic pulsar wind should be very efficient. 

Loading of supersonic flows greatly reduces the strength of the termination shock 
(the change of the Lorentz factor of the flow), but not the kinetic energy flux through it. 
This may affect the efficiency of particle acceleration and PWN luminosity both in radio 
and X-rays. Indeed, Williams et al. (1999) considered the strength of termination shocks 
in mass loaded isothermal flows and found that shocks become weaker. They have also 
suggested that weakening of the shocks may contribute to the radio quietness of some 
wind-blown bubbles. Here we extend this possibility to PWNs. Unfortunately, at this 
point our understanding of relativistic shock acceleration is not good enough to argue 
whether weaker shocks (in terms of a change of a Lorentz factor) are less efficient at particle 
acceleration. 

It is unlikely on numerical grounds that loading of supersonic flow will slow it down 
to the critical point. If it still happens, then, since for a broad variety of the loading 
profiles the critical point of the flow is a focus and for a focus a smooth transition from the 
super- to sub-fastmagnetosonic flow is not possible, a flow would shock at a Mach number 
> 1. For example, Galeev & Khabibarkhanov (1990) have argues that weakly magnetized 
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non-relativistic mass loaded flows shock at a Mach number M = 2. Similar analysis for 
relativistic flows needs to be done as well. 

We have shown that, contrary to naive expectations, loading of the relativistic shocked 
pulsar wind by particles from the filaments initially would result not in a slowing down of 
the flow, but in acceleration. To see the dynamic reason for this strange behavior we first 
note that even if a subsonic flow is weakly magnetized at the source, it becomes strongly 
magnetized as it reaches its terminal velocity. The total energy flux then consist of the 
Poynting flux P = Air^fib 2 and a particle flux ~ (Mo + M)(3 2 /2 (using the nonrelativistic 
expression for simplicity). The Poynting flux, subject to the requirement to transport 
magnetic flux, ^76 = £/2 v / 7rr, is inversely proportional to the velocity: P = S 2 / ' (3. Mass 
loading increases particle energy flux, so that Poynting flux should decrease and velocity 
should increase. 

Acceleration of the wind would result in development of instabilities that would try to 
destroy the reason of the acceleration: the need to transport the magnetic flux. After the 
instabilities have developed and the reconnection destroyed magnetic flux the flow will be 
allowed to decelerate, in line with the idea of Begelman (1998). This deceleration will be 
very quick, (3 ~ r _7//2 , allowing it to match to the boundary of the PWN. 

The other possible application of the present work relates to the structure of the ram 
pressure confined PWNs. The pulsar motion through ISM constantly brings a new neutral 
material in the wind. The typical scale in this case is the stand-off distance for the bow 
shock 



a (39) 

\ 4irc£n m p Vxs 



and a typical time is r^/V^s where V^s is the neutron star velocity through ISM. Combining 
with ( |3"8"D we find that the ratio of the mass loaded induced shock transition to the stand-off 
distance 



(40) 



r b \l-tij (l + <7 ) 

which is typically smaller than unity. Thus, mass loading is extremely important for the 
structure of the ram pressure confined PWNs (see also Bucciantini & Bandiera 2001). 



I would like to thank Mikhail Medvedev, Vicky Kaspi, Elena Amato and Steve Balbus 
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A. Mass loading of relativistic magnetized shocks 



In this appendix we briefly consider mass loading of relativistic shocks. Assume that 
at the shock extra material with density Ap in the rest frame of the shocked plasma has 
been added to the flow. Then the energy, momentum, mass conservation and induction 
equations are 



r - 1 

r 

[hH - 



a 



-p + P + b 2 ^j ^ 2 +p + b 2 /2 



(Al) 



where [] implies a difference between shocked quantities (denoted with a subscript 2) and 
unshocked quantities (denoted below with a subscript 1). Introducing the ratio of the 
two four-velocities AT = /?i7i//3272> the ratio of the three-velocities 7Z = P1/P2, upstream 
magnetization parameter a as a ratio of the magnetic energy density to particle energy 
density h\ — a (pi + j^jPij , and normalizing we find 



pi {a/" 2 f^(/3 2 2 7 2 (i + ^-^/r) 
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2(T-1) 



TZ 



TIT 



(A2) 



The mass loading term is positive definite. 



For a = relation ( |A2| ) reduces to the known jump conditions for the relativistic 
magnetized shock of arbitrary strength. For example for a strong (J\f ^> 1, pi pi) 
magnetized shock this gives 



02 



1 



4(1 + a) 

(this gives @ 2 — T — 1 for a = 0). 



2(r - 1) + Ta ± J(2(T - 1) + Ta) 2 + 8(2 - r>(l + a] 



(A3) 



Relation ( |A2| ) shows that mass loading starts to affect the properties of a strong shock 
when a ~ M 2 p\. A full investigation of the relation ( |A2|) is beyond the scope of this paper. 
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B. Non-relativistic mass loaded flows 



Keeping in mind possible applications of the above results to stellar outflows, below we 
briefly consider non-relativistic mass loaded magnetized winds. The governing equations 







r 2 ( pv 2 /2 + ^—[P + h ' 
\d r [r 2 (pv 2 + h 2 /2)] +3 r p = 



-d r [rhv] = 



-A 



r 2 pv 



R, 



may be simplified if one introduces energy, mass and magnetic fluxes, 
C = v ^h 2 + ^—p + pv 2 /2^j,F = vp, K = vh 
which obey the equations 

d r C ■ 

with solutions 



d r F=-—+R, d r K 
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(B4) 



Introducing a fast magnetosonic wave phase velocity, 



J f 



h 2 Tp 
— + — , 
P P 



(B5) 



and eliminating K, in favor of Vf we get the equation for the evolution of velocity: 

,2 _l Oo,2\ p — 1 



{ V2 - V )) ^ V 



r-l2£- T{v 2 + 2vj) 



-Rv' 



(B6) 



The first term on the rhs represents the effects of pressure on the velocity; it is always 
positive since p = (T - l)/(2r«(2 - r))(2£ - T[v 2 + 2v 2 )) > 0. The negatively defined 
second term is proportional to the rate of mass loading. 

When R = the evolution of the magnetized outflows can be integrated analytically. 
Treating velocity as an independent variable we find from the eq. (|B~6|): 



r ocv 



(2-r)/(2(r- 



(2E v - M v 3 - 2£' 



2 x-l/(2(T-l)) 



(B7) 
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Which shows that a terminal velocity at r — > oo is determined by the third order equation 
for v which involves the energy, mass and magnetic fluxes at the source. 

It is convenient to introduce velocity vq (the terminal velocity of the super- 
fastmagnetosonic flow) and <To (magnetization parameter) by the following relations 

M v 2 M v 2 



E, 



o 



2(1 - ao) 2 



2 = MqvIoq M y3g 

2(1 - ao) ~ 2 1 ' 

When mass loading is not important the fast sound speed expressed in terms of v and 

a is 

2 r - 1 / 2 2\ . 2 -Tvfco . . 

v f = — {y -v) + — — (B9) 

For weak magnetization cr <C 1 the fast sonic flow with v = v f is located at 



(r - i)/(r + i)v = v /2 for r = 5/3. (bio) 

From the eq. ( |B7[ ) it follows that unmagnetized supersonic flow reaches a terminal 
velocity vq (and zero pressure), while subsonic flow slows down as v ~ 1/r 2 — >• (and 
constant pressure). Magnetized flows at r = oo behave somewhat different: at r = oo they 
satisfy 

v 3 - vvl + v^ao = (Bll) 
Eq. (pil|) has two real solutions for (T < 1 corresponding to super- and sub- 



fastmagnetosonic branches. For a small but finite magnetization parameter, the 
super-fastmagnetosonic branch remains almost unchanged, reaching v ~ vq at r = oo, while 
the sub-fastmagnetosonic branch asymptotes to a constant velocity determined by total 
magnetic and energy fluxes: v ~ v Q ao = £ 2 /E (for cr 1). Thus, magnetized flows cannot 
slow down to zero velocity since they have to transport magnetic flux. 

Both sub- and super-fastmagnetosonic branches at large radii have p ~ r~ 2 , p ~ r~ 2T 
and h ~ r _1 , so that the plasma (3 parameter (inverse of magnetization) decreases with 
radius j3 = 2p/h 2 ~ r _2 ^ r_1 \ Both super- and sub-fastmagnetosonic branches expand to 
zero pressure at infinity. Flows which were weakly magnetized at the source (large /3) 
become strongly magnetized as they expand. In both cases magnetized flow have zero 
internal pressure at infinity, so that the terminal fast sound speed is equal the shear Alfven 
wave velocity: 

2 B 2 VqOq ( ^p- at the fast branch v ~ vq , , 

v A = - — = — — ~ < £ (B12J 
4ti"P 2,v I at the slow branch v ~ fo<To 



-20- 



Next we consider special points of the flow. The rhs of the eq. (|B"6|) becomes zero at 



( ^ + \ Rr + I [ Rr*dr) ^ = -1 + % - ^ (B13) 
\2(r - 1) r 2 J J M v 2 V ; 

where parameterization (p8| ) was used. This equation has solution if the maximum on 
the rhs (which is reached at v = 3t>oO"o/2) is larger than 0, which limits magnetization 
parameter to a < 2/(3v / 3) ~ 0.384. 

The condition v — Vf, given by 

(i + r) [i + J ^ j = (r - + (r - 2)-^, (Bi4) 

can be satisfied at a given radius only for one particular value of v. 

To proceed further we assume that mass loading has a power law dependence on radius: 
R~ r - n . Eqns flB13|) and (|B14|) then become 

( 5 + r -n(r+l)) t; 2 vjap 

2(T-1) ^ i; 2 v 



(i + D(i + c) - -(r-i)--(2-ry = o (bis) 

where a ratio of the mass flux due to loading to the initial mass flux 

47r f i?r 2 <ir 47ri?r 3 

C = — —■ = - B16 

M (3-n)M 

has been introduced. Below we assume that £ is an increasing function of radius and that 

n < (5 + r)/(r + i). 

The system ( |B15| ) determines the location ( and velocity at the fast sonic point. It can 
be resolve as an implicit function of v: 

v ((3 - n) v 2 (1 + T) 2 + (-1 + T) (1 + n - (3 - n) Y) v 2 ) 

a ° = ( 3 (3-r)r-n(2-r) (i + r)) v 

c _ 2(-l + r)(-3^ + V) 

^ v 2 (3 (3-T) Y-n (2-T) (1 + T)) 1 ' 

see Fig f|. In the absence of magnetic field, o"o = 0, relations (|B17|) can be resolved explicitly 
for the location and velocity at the fast sonic point: 

4 



c 



T(3 — n) — n — 1 



(r-l)(r(3-n)-n-l) 
\ (l + r) 2 (3-n) V ° (B18) 
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For non-vanishing cr the location of the sound point is pushed to smaller radii (and larger 
v), reaching £ = at a = 2/3\^3 and v = Vo/y/3 independent of T and n. In particular, 
for the simple case of a homogeneous loading, n = 0, and T = 5/3 the fast sonic point is at 
£ = vl - 3v 2 /(5v 2 ), a = 2v/5v (-1 + 8v 2 /v%). 
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Fig. 1. — Loading of supperrelativistic flows. Flows with higher initial Loretz factors and 
with higher magnetization parameters are affected by the mass loading at smaller radii. Flows 
with the same <7o approach the sonic point at weakly relativistic velocities (not shown) along 
the similar trajectories. Flows with higher <r experience a mass loaded shock trasition at 
smaller radii than the flows with smaller <j . 




Fig. 2. — Phase portrait of relativistic magnetized mass loaded flows (n — 0, T = 4/3). 
Left: Unmagnetized case <7q = 0. Right: Magnetized flow with ao — 1. Thick lines are the 
critical solutions which start at 7 = 70 and (5 — a"o/(l + (Jo)- 
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Fig. 3. — Phase portrait of non-relativistic mass loaded flows (n — 0, Y — 5/3). Left: 
Unmagnetized case o"o = 0; focus is located at ( — 1, v — vo/V& (compare with Smith 1996). 
Right: Magnetized flow with a = 0.1; focus is located at ( = —.86, v = 0.44t> . Thick 
lines are the critical solutions which start at v = v o and v = <Jo v o at ( = and, in the case 
o"o — 0, asymptote to zero velocity at infinity. Dashed lines are the fast sonic lines, where 
v = Vf. Dot-dashed line is Alfven line, where v equals Alfven velocity. At the point where 
fast sonic line intersects Alfven line the pressure is 0. Characteristics intersect the fast sonic 
line vertically; no physical solution can extend beyond such point. Physical solutions start 
below the fast critical solution and, in the magnetized case, above the low critical solution. 



-24- 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Fig. 4. — Location of the fast sonic point for non-relativistic flows for different values of cr . 



